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ABSTRACT 

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum 
survey at 33" effective resolution of 170 square degrees of the Galactic Plane 
visible from the northern hemisphere. The BGPS is one of the first large area, 
systematic surveys of the Galactic Plane in the millimeter continuum without 
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pre-selected targets. The survey is contiguous over the range —10.5 < / < 90.5, 
|6| < 0.5. Towards the Cygnus X spiral arm, the coverage was flared to |6| < 1.5 
for 75.5 < / < 87.5. In addition, cross-cuts to |6| < 1.5 were made at / = 3, 
15, 30 and 31. The total area of this section is 133 square degrees. With the 
exception of the increase in latitude, no pre-selection criteria were applied to the 
coverage in this region. In addition to the contiguous region, four targeted regions 
in the outer Galaxy were observed: IC1396 (9 square degrees, 97.5 < / < 100.5, 
2.25 < b < 5.25), a region towards the Perseus Arm (4 square degrees centered on 
/ = 111, b = near NGC7538), W3/4/5 (18 square degrees, 132.5 < / < 138.5) 
and Gem OBI (6 square degrees, 187.5 < / < 193.5). 

The survey has detected approximately 8400 clumps over the entire area to 
a limiting non-uniform 1-a noise level in the range 11 to 53 mJy beam~^ in 
the inner Galaxy. The BGPS source catalog is presented in a companion pa- 



per (Rosolowsky et al. 2010). This paper details the survey observations and 
data reduction methods for the images. We discuss in detail the determination 
of astrometric and flux density calibration uncertainties and compare our re- 
sults to the literature. Data processing algorithms that separate astronomical 
signals from time-variable atmospheric fluctuations in the data time-stream are 
presented. These algorithms reproduce the structure of the astronomical sky over 
a limited range of angular scales and produce artifacts in the vicinity of bright 
sources. Based on simulations, we find that extended emission on scales larger 
than about 5!9 is nearly completely attenuated (> 90%) and the linear scale at 
which the attenuation reaches 50% is 3'.8. Comparison with other millimeter- 
wave data sets implies a possible systematic offset in flux calibration, for which 
no cause has been discovered. 

This presentation serves as a companion and guide to the public data re- 
leas^ through NASA's Infrared Processing and Analysis Center (IPAC) Infrared 
Science Archive (IRSA). New data releases will be provided through IPAC- 
IRSA with any future improvements in the reduction. The BGPS provides a 
complementary long-wavelength spectral band for the ongoing ATLASGAL and 
Herschel-SPIKE surveys, and an important database and context for imminent 
observations with SCUBA-2 and ALMA. 

Subject headings: methods: data analysis - surveys - ISM: clouds - submillime- 
ter: ISM - stars: formation - stars: massive 



http : / / irsa . ipac . caltech . edu/Missions/bolocam . html 
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Introduction 



Millimeter- wavelength continuum surveys of the Galactic plane provide the most efficient 
way to find molecular clumps that are the likely formation sites of massive stars and star 
clusters. The development of detector arrays has made blind surveys of large areas possible. 
Such surveys bypass the need for selection based on the presence of embedded stars or 
star clusters, infrared sources, masers, or radio continuum emission. In particular, they 
can locate molecular clumps before stars form, providing vital information on the initial 
conditions of star formation. These surveys can also provide valuable constraints on the 
physical properties of the clumps, especially masses and mean densities, when combined 
with distance information. 

Galaxy-wide surveys are essential for measuring the impacts of the environment on 
clump properties and star formation activity. Do clump properties vary with Galactocentric 
distance, or with location with respect to spiral arms? Do they depend on the level of nearby 
star formation activity? Answering these questions in our Galaxy will provide the essential 
"ground truth" required for the analysis of distant galaxies where individual clouds and 
clumps are not resolved and only galaxy-wide average quantities can be measured. 



Surveys of the galaxy in low-lying rotational transitions of the CO molecule (Dame 



et al. 2001 Jackson et al. 2006) have identified the locations of large molecular clouds and 



traced the transition from atomic gas in the outer galaxy to molecular gas in the inner 
galaxy. Surveys in CO have detected a strong preference for GMCs to form along spiral 
arms (Stark & Lee 2006). Recent improved measurements of kinematic distances (e.g. Pohl 



et al.||2008 ) and accurate distances from VLB! ( Reid et aL||2009 ) have refined our view of the 
spiral structure of the Milky Way. This new increase in data now allows a more complete 
investigation of the relationship between spiral density waves and star formation. 

Locating the early stages of star formation is quite difficult, and often relies on the 
serendipitous location of a cold, dark cloud which appears as an IRDC. While dark clouds 
were known from optically obscured regions, the "discovery" of IRDCs from the MSX (Mid- 



course Space Experiment) satellite was presented by Egan et al. (1998) and Carey et al. 



(1998). IRDCs have now been extensively cataloged with MSX (Simon et al. 2006a) and 



GLIMPSE 8 /xm data (Peretto & Fuller 2009). In the past decade, compelling evidence has 



emerged which suggests that IRDCs are the pre-cursors to massive stars, and therefore, stellar 



clusters (see, e.g. Rathborne et al. (2006, 2008)). IRDCs, however, require the favorable 
viewing condition of being in front of a bright, mid-IR background for detection. Millimeter 
wave surveys, however, provide an efficient means of identifying cold, dense dust throughout 
the Galaxy independent of the Galactic background, and thus may be seen at large distances. 
Moreover, the millimeter-wave emission is optically thin, allowing the properties of the cloud 
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to be determined. Millimeter wave surveys such as BGPS are essential for understanding 
massive star and cluster formation on a Galactic scale, from the inner to the outer Galaxy. 

Studies of nearby clouds have demonstrated that star formation is far from uniform 
over molecular clouds, but is concentrated in unusually dense regions called clumps or cores. 



which fill a small fraction of the area of most clouds (Lada et al. 1991; Enoch et al. 2007 



Evans et al. 2009). Following Williams et al. (2000) and McKee & Ostriker (2007), we use 



the term 'core' to refer to a very dense region destined to form an individual star or small 
multiple star system and the term 'clump' to refer to a region likely to form a group or cluster 
of stars. While clumps forming massive clusters may be even denser than cores forming low 
mass stars ( Mueller et al.||2002 ), in general, clumps are larger and have a lower mean density 



(McKee & Ostriker 2007) 



To delineate clumps from their parent clouds as traced in CO, line emission from 
molecules which are high-density tracers, such as NH3, H2CO, CS, HCN, etc. are extremely 
useful. Spectral lines also provide excellent diagnostics of line-of-sight motions, temper- 
atures, and densities. However, the interpretation of these gas tracers is complicated by 
variations in tracer abundances caused by freeze-out onto grains, sublimation caused by star 
formation, and subsequent complex chemical processing dvan Dishoeck fc Blake|[l998| . Inter- 
pretation is further complicated by uncertainties in optical depths and excitation conditions 
and the impacts of radiation fields and shocks. These complications make the derivation 
from such molecular line data of column densities, masses, and other physical properties of 
star-forming clumps very difficult. A further limitation of spectral line surveys is that the 
tracers of denser gas tend to be weak, and multi-element receiver arrays have only a modest 
number of elements, making blind surveys very time-consuming. Instruments are beginning 
to be developed to address this hmitation, such as NRAO's K-band Focal Plane Arra}|^and 



JCMT's HARP-B (Buckle et al. 2009). 



In contrast, focal-plane arrays containing hundreds of individual bolometers sensitive 
to millimeter and sub-millimeter (sub-mm) radiation are now available. This advance in 
technology enables blind surveys of the Galactic plane, providing a uniform inventory of 
massive star-forming and starless clumps. Catalogs resulting from these surveys will provide 
the data base for subsequent observations in tracers such as hard-to-excite molecular lines and 
higher resolution studies with instruments like CARMA and ALMA. A significant advantage 
of continuum observations of dust at long wavelengths is the low optical depth of the dust 
( Johnstone fc Bally|2006 ), permitting a straight-forward estimate of the mass of the emitting 
region. 



^https : / / safe . nrao . edu/wiki/bin/view/Kbandf pa/WebHome 
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Submillimeter observations of Galactic sources have been conducted by a variety of 
different groups focusing on small regions. The Submillimeter Common-User Bolometer 



Array (SCUBA; Holland et al. (1999)) on the 15-m James Clerk Maxwell Telescope mapped 
many significant regions in the Galactic plane over its 8 year lifetime 



see 



Di Francesco 



et al. (2008) for a summary), but never completed a contiguous survey of the Galactic plane. 
The Bolocam instrument has been used on the 10.4-m Caltech Submillimeter Observatory 
(CSO|^ to map nearby large molecular clouds at 1.1 mm and identify hundreds of cores 
( [Enoch et aL]|2007D . MAMBO on the 30-m IRAM telescope has been used at 1.2 mm to map 



both small, nearby cores (Kauffmann et al. 2008) and more distant regions of more massive 



star formation (Rathborne et al. 2006; Motte et al. 2007). None of these surveys covered a 



significant fraction of the Galactic plane. 

This situation has begun to change dramatically. Recent wide-area surveys have been 
carried out in the submillimeter bands 250 /xm - 500 /xm by the Balloon-borne Large Aperture 
Submilhmeter Telescope (BLAST) ( |01mi et al.|[2009| |Roy et al.|[20lol [Netterfield eTaL][2009 



Chapin et al. 2008). This work provides an unprecedented view of a wide range of star 



forming environments and demonstrates the richness of submillimeter datasets, creating 
great anticipation for the imminent results from instruments aboard Herschel, particularly 



the SPIRE camera. From the ground, the ATLASGAL survey (SchuUer et al. 2009), using 



the LABOCA instrument at 870 fim with 19'.'2 resolution on the APEX telescope in Chile, 
has covered 95 square degrees in the Galactic Plane, —30 < I < 11.5 and 15 < / < 21 with 
\b\ < 1. The survey will eventually encompass —60 < / < 60 with |6| < 1.5. 

The BGPS is the first millimeter survey of a substantial fraction of the Galactic Plane 
in the northern hemisphere. The BGPS maps in the first quadrant overlap with the VLA 



Galactic Plane Survey (VGPS) in Hi and radio continuum (Stil et al. 2006), the Spitzer- 
GLIMPSE and GLIMPSE-II fields ( [Benjamin et al.|2003p , MIPSGAL ( [Carey et al.|2009| , the 
BU-FCRAO i^CO J = 1 ^ Galactic Ring Survey ([Jackson et al.[[2006D, and ATLASGAL 



( [SchuUer et al.|2009[ ), among others. The BGPS coverage will also overlap with the upcoming 



Herschel HiGAL (Molinari et al. 2010) and SCUBA-2/JCMT Plane Surveys (Di Francesco 



2008 ). The BGPS provides a long-wavelength dust continuum data point, complementing the 



spectral energy distribution provided by other surveys. By covering the both the Galactic 
center and anti-center regions with a single instrument, a consistent comparison over the 
widest possible range of Galactocentric radius is possible. The BGPS gives a view of massive 
star and cluster formation throughout the Galaxy, and provides an important finder chart 
for future high-resolution observations with facilities such as CARMA and ALMA. 



■^The Caltech Submillimeter Observatory is supported by the NSF. 
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This paper presents the imaging data from the Bolocam Galactic Plane Survey (BGPS). 
The source catalog is described in a companion paper ( [Rosolowsky et alT (2010); hereafter 
RIO). The outline of this paper is as follows. Section [2] describes the instrument and the 
observations. The data analysis is described in Sections [3] (astrometry), |4] (the map- making 



algorithm), and |5] (flux calibration). Key discussions are in Section 5.3, which explains the 
effect of the data processing on the accurate recovery of extended structure, and in Section 



5.5[ which compares the BGPS flux with other surveys, including a discussion of a possible 
systematic offset. The final image data products and the public release are described in 
Section |6l We conclude in Section [7] with a brief discussion of the broad features of the 
BGPS in comparison with other surveys. 



2. Observations 

We used Bolocamj^ to survey more than 170 square degrees of the northern Galac- 
tic Plane. Bolocam is the facility 144-element bolometer array camera mounted at the 
Cassegrain focus of the 10.4-m mirror of the CSO on the summit of Mauna Kea. We used 
the filter configuration with a band center of 271.1 GHz (hereafter 1.1 mm) and fractional 
bandwidth I\vjv = 0.17 (46 GHz). The detectors are silicon nitride micromesh absorbers 
with NTD germanium thermistors, operated at a temperature of ~250 mK. The Bolocam 
array field-of-view is 7^5, with individual detectors having nearly Gaussian beams of 31" 
FWHM. The spacing between the individual detector beams on the sky is 38", so the focal 
plane is not instantaneously fully sampled. The Bolocam instrument is described in greater 



detail in Glenn et al. (2003). 



The maps presented here were acquired during six separate observing sessions at the 
CSO over the course of two years between June 2005 and September 2007. The observing 
epochs are given in Table[T] Bolocam observations were typically scheduled when T225 > 0.06. 
Between observing epochs, Bolocam was removed from its mount at the re-imaged Cassegrain 
focus and stored warm. Thus the flux calibration and pointing model were re-computed for 
each epoch to allow for variations in the instrument and optics. The final data products 
were aligned to the well-constrained pointing model created from the Epoch V and VI data 
(Section [3]). We found that the flux calibration did not in fact differ significantly between 



epochs; see Section 5.1 for more details. 



Figure [T] shows the coverage of the BGPS. The survey is contiguous over the range 
— 10.5 < / < 90.5, \b\ < 0.5. Towards the Cygnus X spiral arm, the coverage was flared to 



http : //www . CSO . caltech . edu/bolocam 
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\b\ < 1.5 for 75.5 < / < 87.5. In addition, cross-cuts to |6| < 1.5 were made at / = 3, 15,30 
and 31. The total area of this section is 133 square degrees. With the exception of the 
increase in latitude, no pre-selection criteria were applied to the coverage in this region. In 
addition to the contiguous region, four targeted regions in the outer Galaxy were observed: 
IC1396 (9 square degrees, 97.5 < / < 100.5, 2.25 < b < 5.25), a region towards the Perseus 
Arm (4 square degrees centered on / = 111, b = near NGC7538), W3/4/5 (18 square 
degrees, 132.5 < / < 138.5) and Gem OBI (6 square degrees, 187.5 < / < 193.5). The total 
area of good coverage for the BGPS is 170 square degrees. 

Our basic observing strategy was to raster scan a field by moving the primary mirror 
of the CSO. This was done to modulate the astrophysical signal faster than fluctuations in 
atmospheric opacity. Each field was scanned in a pattern with alternating rasters along lines 
of constant /, followed by a series of rasters along lines of constant b. Fields were observed 
several times with such raster scans to improve signal-to-noise (S/N). Thus each point in the 
final map of a field was observed multiple times with the array moving in different directions 
with respect to the astrophysical emission, helping to separate it from atmospheric emission 
(see Section |4] for further details on this process). 

In Epochs I and II, the fundamental observing block was a region 1° x 1°, covered by 
23 rasters along each of I or b. Starting with Epoch III, the fundamental block was changed 
to 3° X 1° for increased mapping efficiency (i.e., less time spent in turnarounds at the end of 
a raster). Each 3° x 1° block was covered using either 23 rasters along lines of constant b, 
or 67 rasters along lines of constant /. In both cases, the spacing between adjacent rasters 
was 162". The total time for observing the entire block was 39 (48) minutes for rasters 
along constant b (/). The data were electronically sampled at 10 Hz along the scan direction, 
slightly higher than the Nyquist rate for the scan speed of 120" s^^. 

These fundamental observing blocks were stitched together to make 29 large, contiguous 
maps. These images define the "fields" used for comparison. Each "field" includes all 
observations covering a region between 1 and 6 square degrees. 

Because the Bolocam detector beams do not fully overlap on the sky, if the array is 
scanned along detector rows, gaps will be left in the resulting map. To ameliorate this, 
we began using in Epoch II a field rotator to adjust the rotation angle of the array with 
respect to the scan direction so that as the scan proceeded the full extent of the FOV in 
the direction orthogonal to the scan was sampled by at least one detector. This is shown 
in Figure |2} Achieving a uniform sampling along the scan-orthogonal direction is further 
complicated by a number of missing bolometers. A simulation was performed to determine 
the optimal angle to rotate the array with respect to the scan direction to account for both 
the effects of the beam spacing and missing bolometers. During the turnaround following a 
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raster, the field rotator is adjusted to this optimum angle for the subsequent raster. Without 
the field rotator, the coverage shows variations of 100% from pixel to pixel in a single raster 
of a field. With the rotator, this is reduced to ~ 40%. 

3. Astrometry 

3.1. Absolute Reference Sources and Pointing Model 

We constructed an absolute reference system for the BGPS by observing bright quasars 
and blazars near the Galactic Plane with small raster maps. The sources were chosen from 
the SMA Submillimeter Calibrator List[^ since they have positions reliably determined from 
interferometric measurements and are point sources at the scale of the Bolocam beam. The 
distribution of absolute pointing sources over the sky is shown in Figure |3j 

Pointing observations were performed approximately once every two hours over each 
night. The elevation offsets showed a deviation from zero which was empirically well-modeled 
by a quadratic function of altitude. No systematic deviation was observed for the azimuth 
offsets. The derived model is shown in Figure |4j An RMS scatter of ~ 6" for the model 
over the entire observed range was achieved, with somewhat worse scatter in altitude than 
in azimuth. 



3.2. Pixel Positions in the FOV 

In addition to the model for the pointing center, it is necessary to empirically deter- 
mine the actual projected pattern of the array on the sky, and measure the rotation angle 
between the focal plane and the sky coordinate system. This is done by making observations 
which track all detectors across a bright source (e.g. a planet) and making maps from each 
bolometer individually. 

In our first data release, the projected positions of the detectors on the sky are made 
using the nominal positions without correction for optical distortion. This results in a slightly 
reduced effective resolution (33" instead of 31"). The effective beam size was measured from 
fully-sampled beam-mapping observations of planets, and is a good description of the blurring 
of the beam for the portions of a field observed by all the detectors. Within about one FOV 
of the field edges, however, only a subset of the array observes the field, and so the effective 



"http: //smal . sma.hawaii . edu/callist/callist .html 
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beam size varies slightly at these positions, and may be asymmetric. Because this distortion 
is asymmetric, it can also result in a small pointing offset. However, this effect is strongly 
mitigated by the low number of hits at field edges: this data is generally flagged out. The 
residual offset is found to be within the overall pointing model error, but caution should 
be used when measuring source locations at the field edges in the first data release. Future 
releases will have the distortion correction applied, removing this effect and improving the 
effective resolution. 



3.3. Relative Alignment and Mosaicking 

A master pointing field was constructed from observations taken in the epoch with the 
best-constrained pointing model for that field (either Epoch V or VI) for every field observed 
by the BGPS. Relative alignment between the master and other observations of the same 
field was performed by finding the peak of the cross-correlation between images. The cross- 
correlation peak was fit with a Gaussian profile and the difference between the Gaussian 
peak and the image center was used as the pointing offset. The offset may be determined 
accurately to within the error in finding the peak of the Gaussian, typically < 1". To create 
the final image, all observations of a field were merged into a single timestream with these 
pointing offsets applied. 

This method of alignment makes use of all the information available in the maps, and 
avoids the ambiguities inherent in using extracted sources to align fields, as the BGPS sources 
are rarely point-like. It further avoids the slightly larger effective beam and loss of peak flux 
density which would result if the maps of individual observations were combined using the 
pointing model alone, where each individual observation would be co-added with the 6" RMS 
uncertainty of the model. In some fields with few sources, particularly those in the range 
1=65 to 75, there was not enough signal to acquire a pointing offset using cross-correlation. 
In these fields, the pointing model was used directly, and so the above mentioned effects 
affect these fields. 



3.4. Comparison to the SCUBA Legacy Catalog 



All SCUBA 850 and 450 fim data has been re-processed in a uniform manner and made 



publicly available as the SCUBA Legacy Catalog (Di Francesco et al. 2008, hereafter SLC). 
These maps allow a cross-check of the accuracy of the BGPS pointing model. We have 
applied the same cross-correlation procedure used to obtain the relative alignment of BGPS 
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observations to compare the BGPS to the SLC 850 /im images. Di Francesco et aL (2008) 
claim a typical pointing accuracy for SCUBA of ~ 3", with offsets as large as ~ 6" occurring 
occasionally. No further adjustment to the nominal JCMT pointing model was performed 
in producing the SLC maps, and thus we assume these numbers are typical. 

There are a number of differences in the Bolocam and SLC maps which might lead 
to difficulties in comparing their astrometry. In particular, different scan strategies were 
used (SCUBA observations were typically taken in jiggle-map mode, which tends to remove 
extended structure), and of course there is the difference in wavelength and beam size between 
the instruments. Nevertheless, the by-eye morphological comparison between SLC 850 /im 



and BGPS sources is generally excellent, as exemplified in Figure 5(a) A cross-correlation 
technique was used in order to minimize these limitations, since it does not require the 
establishment of a single position for an object. 

The SLC has the advantage of having observations over a large portion of the Plane, 
allowing the pointing to be checked (in spots) for deviations as a function of position. To test 
this, 21 of the 29 BGPS fields were used which overlapped with the SLC. For these fields, we 
find Al = l'!8 ± 1'.'2, Ab = — 0'.'4 ± 0'.'8, nearly consistent with zero, and an RMS dispersion 
about the mean BGPS position of 6'.'2, consistent with the errors derived for the Bolocam 



pointing model. The offset for each field is shown in Figure 5(b) We have not corrected the 



Bolocam pointing based on the SLC or any external data set, but have produced our final 



maps with the internally derived model described in Section 3.1 



4. Mapping Algorithm 

At wavelengths near 1 mm, the emission from the atmosphere dominates any astro- 
physical emission. For an atmospheric zenith optical depth of r = 0.1, the sky brightness 
temperature is around 30 to 40 Kelvin at typical zenith angles. This background loading 



determines the instrument calibration (see Section 5.1), but the primary time-variable signal 
is due to few percent fluctuations in atmospheric opacity. The conversion J from Jy beam~^ 
to Kelvin for Bolocam is 



26 2fc Jybeam ^ 



J= 1 X lO^''— = 58 

where k is Boltzmann's constant and Aes is the effective collecting area of the telescope. 
Thus atmosphere fluctuations are lO's of Jy beam~^, as compared to the typical brightness 
of Galactic sources, which lie in the range of 0.1 - 1 Jy beam~^. The essential signal- 
processing problem to be solved by the mapping algorithm, then, is the estimation of the 
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astrophysical emission in the presence of the fluctuating atmosphere signal, without a priori 
knowledge of either. 

For this work, a custom pipeline was developed to address issues specific to the BGPS. 



It is based on, but significantly improves, the method of Enoch et al. (2006) for generating 



maps and characterizing their properties, and is similar in many respects to the methods used 



by other bolometer cameras, e.g.. Cotton et al. (2009) (MUSTANG/GBT) or Kovacs (2008) 



(SHARC-2/CS0). It also incorporates some of the technical developments of Sayers et al. 



(2009) developed for Bolocam data at 2.1 mm. A key element of the reduction is the iterative 



estimation of the atmospheric fluctuations and the astrophysical signal. The atmospheric 
model is developed from a set of principal components of the bolometer signal time series 
under the assumption that the bulk of the correlated signal is atmospheric. After subtraction 
of the atmosphere model, the astrophysical emission is estimated, then subtracted, and the 
process is repeated. 



4.1. Algorithm 

We assume the raw timestream data d for each bolometer (indexed by i) at discrete 
time points (indexed by t) can be written as 

dit = Sit + au + Cit + Bit + Pit + Sit (2) 

where s is the astrophysical signal, a is the median atmospheric fluctuation seen by all 
detectors, c are terms correlated between bolometers that are distinct from a, e are non- 
random signals due to the instrument itself, p are long timescale fluctuations not modeled 
by a or c, and e is irreducible Gaussian noise due to photon fluctuations and detector noise. 
(Each of these terms is discussed in greater detail in the following sections.) 

In the presence of purely Gaussian noise, the maximum likelihood map m can be ob- 
tained by minimizing the goodness-of-fit statistic 

= {d- AmfW{d - Am) (3) 

Here d is the concatenated data from all bolometers (of dimension A^tod=iiumber time samples 
of all detectors). A provides the mapping between a given sample in the time-ordered 
data and the pixel in the sky map to which it corresponds. A is a matrix of dimension 
Ntod X Npix =number of map pixels. Each element of A is either 1 or 0. Note that we are 
considering the data d and the map m as vectors (the mapping to a two-dimensional map 
being implicit in the matrix A). W is an Ntod x Ntod matrix which is the inverse of the 
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covariance matrix of the time domain noise e, 

Nw = {stst'), (4) 
where () denotes an ensemble average over many reahzations of the random timestream noise 



process St- Further detail about the way in which W was estimated is given in Section 4.6 
The result of the minimization is 

m = {A^WA)A^Wd (5) 

Thus, the mapping from data to map is a linear operator, M = {A'^WA)A^W . For com- 
pactness. Equation [5] will be written as M[d] = m. Note that the mapping operator is not 
invertible. However, given a map m and an observing matrix A, there exists a linear opera- 
tion which makes predictions about the observed timestreams, namely d = Am; this will be 
denoted T[m] = d. 

Because the model Am in Equation [3] only includes the celestial contribution to the 
observed data d, direct use of Equation |5] will produce a map containing artifacts due to the 
unmodeled components a, c, e, and p. One approach would be to include templates for these 
terms as additional rows in A and proceed with a simultaneous fit. In general, however, the 
correct forms to use are not known a priori. The goal, then, is to produce a time series for 
each bolometer which as closely as possible approximates s + e so that we may produce the 
best estimate of the astrophysical signal m = M[s + e] = S + N. We proceed iteratively, 
estimating each term in Equation |2] in the order of its relative strength. In the following, we 
denote the best estimate of a time series x in the n*'^ iteration as x*-"-*. 

In the first iteration, we estimate the largest signal, a, according to the method in 
Section 4.2, and assume sfl^ = cf^^ = pfl^ = 0. We then form 

= dit - aj"^ - c[f - 4"^ - A? - Sit + Sit (6) 

This is the best estimate of signal plus irreducible noise at iteration step n, and is made into 
a map 

M[dr^] = (7) 

The current best map m*^") is then deconvolved to provide a relatively low-noise, smooth 
map from which to generate a timestream 

T[P[m(")]] = 4") (8) 

where T> represents the deconvolution operation. At this stage, we can subtract both a^-"'' 
and 4^ from the original data, and estimate c-"\ and p^^\ The iterative process begins 
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again with Equation |6} An example time series showing the successive removal of these 
models is shown in Figure [6] We discuss the estimation of each of the various terms a, c, e, 
and p in the following sections. 

It is useful at each step to produce a "residual map" 

E = M\du-S^-5!i^-i^-i^-f^] (9) 

This residual map serves as a visualization of the progress of the iterative method, since in 
a perfect process it would be map of e. 



4.2. Atmosphere Fluctuation Model 

Since atmosphere fluctuations are the largest signal in the raw data, they are modeled 
first. By construction, the Bolocam beams as they pass through the dominant layer of atmo- 
spheric water vapor a few km above the telescope are still highly overlapped, and therefore 
sample a nearly identical region of atmosphere. The simplest model for the atmosphere 
fluctuations makes use of this fact by assuming that the largest correlated signals between 
detectors are due to to the atmosphere fluctuations; this is the term a in Equation |2} We 
construct a model of the atmosphere fluctuations which are seen in common by all detectors, 
which we denote the "median atmosphere model". In each iteration n, we construct the 
non-atmosphere-fluctuation part of the time series as 

5i^=du-S^-^-^^-p^ (10) 

The median atmosphere model is then simply 

a|"^ = mediani[a-"^] (11) 

Each bolometer's timestream is then fit to the median atmosphere model, i.e., we find r^'^'' 
such that a[^^ most closely approximates r^'^^af^\ This defines the relative gains of the 
detectors as 

E~ (n) ~ (n) 



The Ti converge fairly rapidly. (See the discussion in Section 5.1 



The simple common mode model of Equation 11 does not remove all of the signal 



correlated at zero time lag between the bolometers. Some of this remaining correlated signal 



may also be atmosphere. Sayers et al. (2010) found for the 2.1 mm Bolocam data that 



the subtraction of a second order spatial polynomial across the array accounted for nearly 
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all of their residual noise; however, the 1.1 mm data seem to have additional correlated 
components that are not completely described by such a model. Therefore, to further remove 
the atmosphere (and any other correlated instrument noise) without a detailed physical 
model, we use a principal component analysis (PCA) as given in Laurent et al. (2005). 
While this method was developed for point source detection in deep, high-redshift surveys, 
the technique works equally well for this survey, with the caveat that it removes large-scale 
structure on the size of the array FOV. This effect is mitigated by the iterative process (the 
Laurent et al. (2005) analysis did not iterate). Figure [T] shows an example of the modes 



removed by the PCA model. These are a graphical visualization of the terms c in Equation 

m 



4.3. Instrument Error Signals 

Most of the instrument error signals have characteristic features which aid in their 
removal. In particular, unlike the atmosphere emission and correlated signals described 
above, these signals are not correlated between detector channels. The error signals include 
the following: 



1. Pickup from the 60 Hz AC power. This appears as narrow lines in the power spectral 
density (PSD) of the data. The second harmonic of 60 Hz is aliased via beating against 
the 130 Hz bolometer bias frequency into 10 Hz with sidebands split at ~ 1 Hz. It is 
removed by replacing components at these frequencies in the Fourier transform with 
Gaussian noise matched to the local mean amplitude. 

2. Spikes in voltage due to cosmic ray strikes on the bolometers ("glitches"). An example 
of a glitch in the time series is shown in Figure [8| along with the extracted distribution. 

3. Microphonic pickup due to vibrations of the receiver. The microphonic effect is due to 
a change in the capacitive coupling of the readout wires to the circuit ground, which 
is converted by the high impedance of the detectors into a measurable signal, much 
like a condenser microphone. The most noticeable microphonics occur at the end of 
each scan when the telescope is turning around and the field rotator is adjusting. This 
leads to broad spikes in the time series, whose long decay must be removed from the 
data, particularly during the beginnings of scans. 



Fortunately, most of the AC powerhne pick-up occurs at frequencies where there is 
no astrophysical signal, given the beam size and scan speed (c.f. Figure 1 of Sayers et al. 



(2009)). Thus this error is dealt with by first notch filtering at the line frequencies and then 
low-pass filtering the data. Because of the low correlation with astrophysical signal, this step 
is performed only once and is not iterated. 



Both glitches and microphonic pickup from the scan turnarounds are degenerate with 
astrophysical signal and must be estimated as part of the iterative process. Glitches are 
identified as large excursions from the RMS level after subtraction of the best atmosphere 
and bright source model, and the data there is excluded in subsequent maps. The turnaround 
microphonics are explicitly modeled as decaying exponentials at the beginnings and ends of 
scans. Residuals from this model may also be removed by the PGA cleaning. 

After iteratively subtracting the astrophysical model, each scan also has a fifth order 
polynomial fit in time removed to deal with the longest time scale modes (p in Equation [2]) . 
These modes are not correlated between detectors. They correspond to spatial scales much 
larger than those to which Bolocam remains sensitive after subtraction of the PGA model, 
so their subtraction produces maps with fewer "stripes" , without appreciably degrading the 
signal. 



4.4. Data Flagging 

Due to the large volume of data generated by the survey, it was necessary to develop 
new tools to quickly visualize the data and ensure data quality. A common tool in radio 
astronomy is the so-called "waterfall" plot, which is an image with frequency on the x-axis, 
time on the y, and an intensity proportional to the interferometer visibility amplitude or 
phase. We have used a variant on this to visualize the Bolocam data by displaying bolometer 
number (related to the bolometer position in the focal plane) on the x-axis and time on the 
y, with the intensity given by the bolometer's response (in Jy) at that time. An example 
of this is shown in Figure |9] Anomalies such as the glitch shown can be detected manually 
(with practice) and interactively flagged out. 

An automated flagger was also created that flags out outlier data on a per-map-pixel 
basis. In order to make a robust measurement of the variance of the fluxes assigned to 
each pixel, we used the median average deviation (MAD) over the data in the pixel and 
rejected high and low outliers at the 3-sigma level. Pixels with too little data to compute a 
deviation, i.e. those with < 3 data points, were also flagged out - these scan-edge pixels are 
the dominant contribution to the total number of flagged data points. The fraction of the 
data flagged by all methods was about 0.08%. 
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4.5. Creation of the Astrophysical Model 

The timestream data is made into a spatial map using the pointing data corresponding 
to each time point for each bolometer. The data is weighted by inverse variance across a 
single scan and then drizzled into a map with 7'.'2 pixels using a nearest-neighbor algorithm. 
The nearest-neighbor matching allows the map to be returned to a timestream in the same 
manner, but with the S/N improved by averaging over all hits on a given pixel. It has 
the disadvantage that it accentuates the unevenness of the cross-scan sampling, since each 
timestream point's value is assigned to an area much smaller than the beam. In principle 
this can be addressed by various gridding algorithms, and may be a part of future releases 
of the BGPS data. 



The resulting map is then subjected to a maximum-entropy clean algorithm (HoUis et al. 



1992) using a specified kernel to produce a deconvolved image which is mapped back into 
the timestream to be subtracted (Equation [s]). This has several advantages over using the 
original map directly, including rejecting artifacts smaller than the kernel and decreasing the 
noise of the resulting timestream. It is also better than a "threshold" method, which tends 
to produce discontinuities and include noise outliers. 

In some fields, fiux from some sources spread out over the course of the iterative process 
using certain kernel sizes. The cause of this artifact is not yet well-understood, but using 
a different kernel size was an effective workaround and did little to change the properties 
of the final map. The standard kernel size was a 14'.'4 FWHM Gaussian; a 21'.'6 kernel was 
used when artifacts were encountered. Simulations indicate that this change in kernel has a 



negligible effect on the flux of the source; see Section 5.1 



We also investigated whether the choice of deconvolution kernel produced a noticeable 
effect on the recovered flux. Figure [lO] demonstrates the effects of using various kernel sizes. 
The best results, as determined by the flux recovered in simulations and examination of the 
residual maps, were produced by a kernel smaller than the beam but larger than the pixel 
size used. 



4.6. Noise Estimation and Residuals 

The weighting of the time stream data used to produce a map is given in Equation |4j 
If the noise were uncorrelated between timestream points, then Nu' would be diagonal. If 
the noise is correlated between time samples, but is stationary, i.e., the noise properties are 
time-translation invariant 

Nu'=N(\t-t'\), (13) 
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and the weighting can be done using the Fourier transform of the correlation, which is 
diagonah 

N^^' = {\e{io)\y{io~io'), (14) 

where is the power spectrum of the noise et. Finally, if the noise is Gaussian, then 

the covariance matrix provides all the necessary statistical information about the noise. 

In practice, few of these assumptions hold for the BGPS data: the noise is not sta- 
tionary, due to atmosphere variations, and the unsubtracted astrophysical signal causes the 
residuals from which we estimate the noise to be non-Gaussian. Further, the true remaining 
correlations between detectors due to unsubtracted atmosphere are difficult to estimate. Our 
approach, then, is to use a non-optimal but reasonable weighting, estimating the noise from 
the integral over the power spectrum of the residual timestream for each bolometer and each 
raster separately. This accounts for the non-stationarity but ignores the correlations, which 
are small after the PGA subtraction. 



The residual maps E (Equation [9]) form the basis for the noise estimation in the map 
domain, and also provide a way of estimating the systematic error resulting from imperfect 
subtraction of bright sources. These maps are produced for each region from the data 
residuals after removing the terms s, a, c, p, and e in Equation |2} An example of a residual 



map is shown in Figure 10 In an ideal case, the residual map E represents a map of a 
realization of the underlying irreducible noise e. However, the imperfect estimation of the 
signal results in "ghosts" of bright sources. These features are a guide to the number of 
iterations necessary and an estimate of the remaining systematic error. In constructing the 
BGPS catalog, RIO found that smoothed versions of the residual maps provide a reliable 
means of determining the pixel-to-pixel error which accounts for the local variations in the 
noise, including those due to such artifacts (c.f. RIO Figure 2). Because observing conditions 
varied widely during the survey, the RMS noise level varies over the BGPS. However, the 
noise within a given 1° x 1° field is fairly uniform in the absence of very bright sources. 
The noise level varies between 11 and 53 mJy beam~^ in the inner Galaxy. We show the 



variation of the depth as a function of Galactic longitude in Figure 11 This noise level 
interacts with the source density to produce the completeness of source extraction. The 
variation of completeness as a function of Galactic longitude, determined by simulations of 
the source extraction, is shown in RIO Figure 9. 
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5. Calibration and Photometry 

5.1. Flux Density Calibration 

The absolute flux calibration is derived from observations of Mars and Uranus (the 
"primary calibrators"). The millimeter- wave flux of these planets is known to ~ 5% ( [Orton 



et al. 1986 Griffin & Orton 1993). The fluxes from these models were extracted using the 
JCMT's FLUXE^ program with a central frequency 271.1 GHz and bandwidth 46.0 GHz 
(Table 11. 



Calibrations at (sub)millimeter wavelengths are strongly affected by atmospheric opac- 
ity corrections. Further, the detector responsivity of Bolocam's bolometers is a non-linear 
function of the mean atmospheric loading. To address both of these problems and relate 
observations of the primary calibrators to observations of the BGPS fields, we make use of 
the following relation. The calibration C, referenced to the detectors, is given by 



C 



Jy. 



7^(r)r7Aexp(-r)A^/ (15) 



where TZ is the bolometer responsivity ([V/W]), rj is the system optical efficiency, A is the 
effective telescope collecting area, Au is the bandwidth, and r is the line-of-sight, in-band 
atmosphere opacity. Under the assumption that the only power variation on the detectors 
is due to the power from the atmosphere (i.e. that astronomical sources are faint relative to 
the atmosphere), which may be parametrized by r, C is a single- valued function of r. We 
have used the measured potential difference across the detector thermistor Vdc as a proxy 
for r, since the bolometer resistance is a single-valued function of loading. This quantity is 
monitored continuously for all observations. Note that this calibration curve folds in both 
the effects of changing atmospheric transmission and the changes in the detector response 
with optical loading. 

The calibrator field observations were obtained and reduced differently from the science 
fields. Calibrators were observed in ~ 8' radius fields with two sets of observations, scanned in 
orthogonal directions. In these small fields, the iterative process tended to diverge with more 
iterations because of instabilities in the deconvolution procedure. Therefore, the calibration 
curve was determined from non-iteratively-mapped calibrator maps. A future work will 
explore applying a more consistent application of the processing to both calibrator and 
science fields. 

A fit to the observed values {Vdc vs. C) for the primary calibrators was performed 



http : //docs . j ach . hawaii . ed.u/star/sun213 . htx/ sun213 . html 
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for each epoch separately. The agreement between epochs was good, so a single combined 
calibration was used for all data. The calibration curve is shown in Figure 12 The resulting 
error on the calibration curve fit is less than 8% (statistical) over the observed range of r. 

The above flux density calibration only accounts for the average calibration of bolometer 
Volts to Janskys. It is also necessary to account for the variation of bolometer response across 
the focal plane (often referred to as "flat-fielding" ) . We do this by monitoring the response 
to the atmosphere emission in all bolometers (astrophysical signal is negligible in comparison 
to the sky except towards SgrB2). Being in the near field, the atmosphere is common to all 
detectors, and thus this serves as a common relative reference to calibrate out variation in 
the individual detector responsivities and optical efficiencies. (For more detail regarding the 
properties of atmospheric noise above Mauna Kea as inferred from the Bolocam data, see 



Sayers et al. (2010).) The change in relative response with loading for a typical detector is 



shown in Figure \T2\ The actual application of the method to the time series of the bolometer 
data is shown in Figure 13, where it is clear that correlation of the timestreams is indeed 
improved by the application of a single multiplicative factor (the of Equation 12). 



5.2. Point Spread Function, Aperture Corrections, and Surface Brightness 

Calibration 

The Bolocam point spread function (PSF, or "beam") is measured using the planets 
Uranus and Neptune, which are nearly point sources for Bolocam. By stacking all obser- 



vations of planets, we obtain a high S/N profile of the beam; this is shown in Figure 14 
The main lobe is fit to a Gaussian profile. This yields the effective beam FWHM of 33", 
corresponding to a sigma of 14'.'2 and equivalent tophat radius of 19'.'8. The corresponding 
solid angle is 2.9 x 10~^ steradians. 

The calibration from Jy beam~^ into surface brightness (MJy sr~^) is obtained once the 
beam area is determined. For the BGPS, the conversion factor is 

1 Jy beam"^ = 34.5 MJy ster"^ = 0.042 Jy pix"^ (16) 

for the standard 7'.'2 pixel size. The uncertainty in the mean beam diameter is 3%, leading to 
a corresponding 6% uncertainty in the beam area (and in the surface brightness calibration). 
This factor is in addition to the uncertainty in the best fit Voc versus C curve. 

Due to the sidelobes of the beam, aperture photometry may not capture all of the flux 
from a point source if the aperture is too small. RIO used a standard set of apertures, and 
the multiplicative correction factors to be applied to fluxes obtained from these apertures for 
a point source are given in Table |2} The correction is only significant for apertures < 80", 
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beyond which effectively all of the flux is included. For extended objects, the situation is 
more complicated because of the effects of the mapping algorithm. This is discussed further 
in the next section. 



5.3. Effects of the Data Reduction Procedures on Photometry 

To do accurate and meaningful photometry on the Bolocam maps requires a good esti- 
mate of the noise (Section |4.6 ) but also an understanding of the spatial filtering imposed by 



the observing strategy and the cleaning and mapping of the data. The additional complexi- 
ties of defining "sources" for the complicated morphologies actually observed are discussed 
in Rosolowsky et al. (2010). Here we restrict our discussion to the limitations on performing 



photometry imposed by the data reduction. 

The fundamental feature of the data is that there is a degeneracy between large spatial 
scale astrophysical emission and the atmosphere and other signals with long spatial or tem- 
poral variations. Consequently, there are limits to the ability of the algorithm to separate 
the signal s from the atmosphere a and time-correlated terms c in Equation [6j Inevitably, 
some of the signal which should be present in s is mixed into these terms and is thus not 
present in the final map m. Because s is generated from a deconvolution algorithm, it has 
non-zero mean by construction, but importantly, it does not contain an unbiased estimate of 
the astrophysical signal on all spatial scales: large spatial scales are preferentially attenuated. 

The use of a constant factor to convert to surface brightness implicitly assumes that the 
response to all angular scales larger than the beam is the same. In fact, the data reduction 
process acts differently on sources of different angular extent. In performing photometry, 
the flux density in an aperture should in fact be written as 

F = fi^src) [ s/n,dn (17) 

J aper 

where S is the flux density in each pixel in Jy beam~^, Qj, is the beam solid angle in steradians, 
and f{flsrc) is a correction factor depending on the true angular size of the source. The factor 
/ can deviate signiflcantly from unity for extended sources, and as Qgrc is generally unknown, 
this leads to the largest systematic uncertainty in obtaining photometric fluxes. 

To determine the factor f{Qsrc) (the flux recovered as a function of angular scale of 
the source), we performed a series of simulations. Gaussian sources with a range of peak 
flux densities and sizes were inserted on top of a background consisting of real BGPS time 
series data with the astrophysical source model removed. The simulated time series were 
then processed identically to the real data. The recovered fluxes were measured in ellipses 
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that included single well-separated sources. Figure 15 shows the fraction of flux recovered 
in the map as a function of source size for a range of different PCA components subtracted 
in the cleaning. The final reduction used 13 PCA components because this produced a rea- 
sonable compromise between attenuating extended structure and cleaning the atmospheric 
contribution. Note that this is a substantially more aggressive cleaning than that used in 



the reduction of the Bolocam observation of the Spitzer Cores-to-Disks fields (Enoch et al. 



2006 Young et al. 2006; Enoch et al. 2007), largely because our mean integration time per 



pixel was ^seconds, whereas theirs was ~minutes. The recovered flux does depend on the 
source size, but importantly does not depend on the source strength except for the strongest 
sources (e.g., SgrB2). 

While the precise effect of the flux reduction requires a full simulation of the data 
reduction pipeline to ascertain, the flux reduction effect can be reproduced with a relatively 
simple prescription which is intuitively motivated. Because the atmosphere model a and 
PCA model c remove from the timestream "signals" which are correlated spatially across 
the array at the same time, this effectively acts as a spatial filter on the map, removing 
spatial scales comparable to the array FOV. We have found that a simple "brick wall" filter 
in Fourier space which nulls all modes below 0.1 inverse arcminutes (spatial scales larger than 
10') reproduces very well the attenuation seen in the simulation. The comparison between 
the flux attention produced by this filter and that found from simulation is shown in Figure 
15 The implication is that the BGPS is not sensitive at all to scales larger than 10'. 



5.4. Bandpass Effects: Color Corrections and Line Contamination 

The width of the Bolocam passband is such that the effective band center, or equiva- 
lently, the flux density referred to a fixed band center, will change somewhat depending on 
the source spectral energy distribution (SED). Atmospheric transmission variations may also 



affect the effective passband. Figure 16 shows the passband. Bolocam is largely insensitive 
to variations in the strong water absorption lines at 183 and 325 GHz. We compute color 
corrections due to changes in the source spectrum in Appendix |A| 

Note that we have made no correction for the contamination of the continuum flux 
densities by emission from lines. The Bolocam passband was specifically constructed to 
exclude the strong ^^C0(2 — )■ 1) emission line at 230 GHz. We can estimate the fractional 
contribution of a line to the continuum as 



L _ J J T(j){u)t{u) du ^ JTAvUc 
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where T is the hne strength in Kelvin, 0(z/) is the hne shape (of width Af in km s~^), t(z/) 
is the Bolocam passband (with equivalent bandwidth Az/ and bandcenter z/^), and J is given 
by Equation [ij The approximation holds for lines in the center of the band. The scaling 
values are chosen to be typical of the line-to-continuum actually seen. Thus any given line is 
likely to at most contribute a few percent to the measured flux, but in the most extreme star 
forming environments (e.g., Sgr B2) we may expect that lines may contribute a substantial 
fraction of the total measured flux density due to the integrated effect of the lines. For an 



example, see the lower panel of Figure 16 



For completeness, we note that towards Hll regions, free-free emission may represent a 
non-trivial fraction of the millimeter flux density in addition to the thermal dust continuum. 
This effect is of course best assessed on a case-by-case basis by comparing the Bolocam maps 



with radio surveys, e.g., the VGPS (Stil et al. 2006). 



5.5. Comparison to Other Surveys 



Comparison of the Bolocam flux densities is easiest for an instrument with the same or 
similar bandpass, since the additional complication of the source spectral index is avoided. 
For this reason, we choose to make an initial comparison of the accuracy of our calibration 
with the MAMBO and SIMBA instruments, since their passbands at 1.2 mm are quite close 
to that of Bolocam. MAMBO operates on the IRAM 30-meter diameter telescope on Pico 
Veleta in Spain, and SIMBA operated on the 15 m SEST antenna prior to that facility's 
closure. An important future cross-check will be comparison of the BGPS flux density 
calibration with SCUBA 850 /xm or ATLASGAL 870 /im results. 



The study of the Cygnus-X region by Motte et al. (2007) (hereafter M07) presents one of 
the largest fields in the Galactic plane that has been mapped with MAMBO which overlaps 
with the Bolocam survey fields. M07 used both the 37 element MAMBO and 117 channel 
MAMBO-2 instruments. The FWHM beam size for the observations was 11". The MAMBO 
beam model resulting from the chopping-scanning observing strategy is described in the 



of 


Motte & Andre 


( 


2001) 


Kauffmann et al. 


( 


2008 



calibration, with an estimated uncertainty in the fiux of 20%. The effective passband of the 
images presented by M07 is about 90 GHz centered at 260 GHz and thus includes the bright 
J = 2-1 CO transition (see Figure 16). 
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A survey of 2 square degrees towards ^ = 44° was conducted by Matthews et al. (2009) 
(hereafter M09). These images have 22"resolution, and unhke the M07 survey, which used 
chopping and scanning simultaneously, the SIMBA data were obtained in fast scan mode 
only. 

To compare the BGPS images to the M07 and M09 images, we first obtained the data in 
FITS image format. The M07 data is publicly availably The M09 data was kindly provided 
by H. Kirk (private communication). All data were re-gridded to a common projection and 
coordinate system grid with 4" pixels using the program Montag^ We then converted to 
MJy sr~^ in a given pixel based on the beam and pixel sizes. Finally, the Bolocam map 
was convolved with either an 11" (M07) or 22"(M09) FWHM Gaussian and the M07 and 
M09 maps are convolved with a 33" FWHM Gaussian so that the effective resolution of the 
two pairs (BGPS-M07, BGPS-M09) is the same. The convolution preserves the total flux in 
the map, but changes the aperture in which that flux appears. This results in maps which, 
were the calibration and observed features in each map identical, can be compared on a 
pixel-by-pixel basis. 



The plots of pixel values versus each other are shown in Figure [T7j The BGPS values are 
indeed correlated with both MO 7 and M09 over a range of about two orders of magnitude, 
but with a large amount of scatter. We attribute the scatter to both noise but primarily to 
the difference in the handling of extended structure between the different surveys. However, 
the average value of the correlation coefficient deviates significantly from 1 in comparing to 
both M07 and M09, with a multiplicative factor of 1.21 to 1.51 required to make the BGPS 
fluxes match the other surveys (or 0.66 to 0.83 to make the other surveys agree with the 
BGPS). This factor also depends on the range over which the fit is done, as given in Table [3] 
Including all pixels down to 3 MJy sr~^ produces systematically larger required correction 
factor. Using the average value (to account for variation of the factor with brightness) yields 
a factor of 1.35. For a dust spectral index of ~ 3.5, the 1.2 mm measurements should 
have produced a flux that is about 12% lower than that measured by BGPS (see Table |4]). 
Assuming this is a typical spectral index for these clumps, we therefore note that one 
must multiply by a factor of 1.35 x 1.12 = 1.5 (with a range of plus or minus 
0.15) the Version 1.0 images and catalog released through IPAC to obtain good 
agreement with these surveys. 



' ftp : //cdsarc . u-strasbg . f r/pub/cats/ J/A+A/476/1243/f its/ 

^This research made use of Montage, funded by the National Aeronautics and Space Administration's 
Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Num- 
ber NCC5-626 between NASA and the California Institute of Technology. Montage is maintained by the 
NASA/IPAC Infrared Science Archive. 
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The source of this discrepancy remains unclear. It is unhkely to be due to hne contam- 
ination in the MAMBO / SIMBA bands. Line contamination could increase the continuum 



flux by at most a few percent for MAMBO (see Equation 18), while being negligible for 
Bolocam. For most sources, it is also not obvious that the discrepancy is due to filtering 
of extended structure. The filtering does not reach the requisite 50% needed to explain the 



difference until the source FWHM reaches 230" (Figure 15). However, the simulations were 
run on Gaussian sources, and it is still not certain whether the filter function differs for other 
source shapes. The effective filter functions of the SIMBA and MAMBO surveys are also 
not well understood. 

Work is ongoing to thoroughly investigate all potential sources of this discrepancy with 
further more realistic simulations, a new treatment of atmosphere removal and application 
of distortion corrections. In addition, the processing for bright flux density calibrators in the 
current pipeline is different from the processing of the science fields, and attempts are being 
made to process the different observation types with an identical pipeline. This analysis will 
be discussed in a future work. 



6. Final Maps and Data Release 

The final maps are produced by co-adding all observations of a given region. The maps 
were processed in 29 separate pieces. The maps are made in Galactic coordinates using a 
plate carre (FITS header CAR) projection. (This is the same projection used by the Spitzer- 
GLIMPSE and BU-FCRAO GRS surveys.) By placing them all on the same equatorial grid, 
mosaicking the maps together is straightforward. As these maps are near the coordinate 
system equator, the difference between truly equiareal pixels and the pixels used is at most 
0.4%, even 5° out of the plane (as for IC1396). The pixel size is 7'.'2, chosen to be much 
smaller than the Bolocam beam. The maps are written to standard FITS files. An example 
of the FITS header for the BGPS images is given in Appendix [B| 

All processed maps are available through IPAC|^ IPAC provides a cutout service for the 



images as well as a searchable version of the catalog provided in Rosolowsky et al. (2010). 



The Version 1.0 release includes the following images, each covering the same regions: 

1. Calibrated maps for all regions observed (MAP). The maps are standard astrophysical 
maps subject to the caveats previously mentioned, the most important being the spatial 



filtering function (Section 4.6). 



http : //irsa. ipac . caltech. edu/Missions/bolocam.html 
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2. Maps of integration time per pixel (NHITSMAP). The nhitsmap shows how many 
recorded data points from the timestream have been assigned to each pixel; it is pro- 
portional to the dwell time per pixel (each 'hit' is 0.1s). 

3. Maps of the residual time series (NOISEMAP). The residual map is the map of the 
model timestream subtracted from the data timestream. 



Examples of each are given in Figure 18 Future releases may contain additional data 
products. 



7. Discussion 



The BGPS has observed a 170 square degree area at an effective resolution of a Gaussian 
PSF of 33" FWHM (corresponding to an area of 2.9 x 10^^ ster). The positional uncertainty 
of the maps is 6" RMS. The maps are calibrated into Jy beam~^ using Mars as a primary 
flux standard, and accounting for atmosphere opacity variations in real time. The BGPS is 
well suited for studying source structure smaller than 3 arcminutes in size. The processing 
of the maps removes more than 90% of flux for features with extents larger than 5'.9, and 
attenuates the aperture flux of structures extending to 3'8 by 50%. We have compared our 
flux densities to those of M07 and M09 and find that we need to multiply our flux densities 
up by a factor of 1.5 to match. We have investigated possible sources for the discrepancy, 
but do not understand the source of the difference. For now, we recommend multiplying the 
flux densities in the maps presented here and in the RIO catalog by the factor 1.5 ± 0.15 
to obtain consistency with other data sets. We are continuing to explore the source of 
this discrepancy, particularly in the calibration and filtering stages of the data pipeline. The 
latest data releases and updates on the discrepancy can be found on the BGPS IPAC website. 



Maps from the entire BGPS are presented in Figure 19 Figures [5] and 19d show 



the generally excellent morphological agreement between BGPS maps and those at other 
(sub) millimeter wavelengths. Filamentary structures with aspect ratios up to the maximum 
spatial dynamic range of the survey (354"/33" ~ 10) are present throughout the Galactic 
plane. The most crowded fields consist of frothy and clumpy structure. The majority of 
sources are at least moderately resolved. We suspect that higher resolution observations 
will resolve BGPS clumps into clusters of protostellar objects, as has been seen in Serpens 



(Enoch et al. 2008; Testi & Sargent 1998), OMC 1 (Beuther et al. 2004 Johnstone & Bally 



1999), and S255N (Cyganowski et al. 2007) where bolometer observations at low resolution 



have been complemented by interferometer observations at high resolution. 
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Figure 20 shows a comparison of the BGPS with Spitzer-GLIMPSE 8 /im (Benjamin 



et aL||2003D and VGPS 20 cm continuum images ( |Stil et aL||2006D . In Figure |20| IRDCs 
appear as the bright white clumps, often filamentary, near the mid-plane. BGPS sources 
are clearly associated with some of the Hll regions detected by the VGPS; a zoom in of this 
region is shown in Figure |2Tj The comparison with 8 /im IRDCs is particularly striking, 
with a high degree of correspondence in their morphologies. Detection of an IRDC requires 
that the cold, dense cloud be on the near-side of a bright mid-IR background, whereas BGPS 
sources are subject only to a sensitivity limit. Comparison of the GLIMPSE and BGPS data 



(see Figure 22) show that some BGPS sources are clearly associated with IRDCs, a few are 
associated with IR-bright features (likely due to the association of free-free emission in the 
millimeter and PAH emission at 8 /im), and some have no clear IR-counterpart. The BGPS 
sources which are neither bright nor dark at 8 /xm are probable candidates for "IRDC-like" 
clumps (cold, dense, potential proto-clusters) on the far-side or outer edges of the Galaxy. 
Thus, the BGPS could provide candidate pre-clusters independent of the favorable viewing 
angle of being on the near-side of bright mid-IR emission. An exploration of this population, 
as well as a detailed analysis of the physical properties and star formation tracers of IRDCs 



corresponding to BGPS sources may be found in Battersby et al. (2010). It is worth noting 
that many features in the BGPS maps have no corresponding counterpart in the Spitzer or 
radio image. 

There are clear maxima in the emission filling factor in the Galactic Center and in the 
region from / = 23 to / = 31, the line of sight through the densest part of the molecular 
ring. Additional areas of high source density are seen at / = 10 and / = 13, as shown in 



Figure 23 This result is in contrast to CO filling factors which are ~ 1 for / < 40 (Dame 



Figure 


23 




et al. 


2001 



et al. 2001; Jackson et al. 2006). In parts of the plane where only one CO emission line 



is seen along any given line of sight (/ > 40), the integrated emission is much lower, but 



there are still detections of CO in every pixel of the Dame et al. (2001) map in the latitude 
range —0.5 < b < 0.5. This behavior is consistent with the interpretation that the BGPS 
preferentially detects denser material than CO. The connection of the millimeter emission 
with star forming material awaits the detailed comparison to other star formation tracers. 
BGPS clumps are, however, associated with well known star forming regions: the brightest 
clumps visible in our maps are Sgr B2, G34.3-h0.15, W 51, W 43, W 49, and M 17. These 
associations suggest that other bright sources of millimeter emission will prove to be massive 
star- or cluster-forming regions. 

In some cases, the distance to BGPS clumps can be determined by matching Galactic 
Ring Survey ^'^CO and 1.1 mm continuum morphology (Simon et al. 2006b Jackson et al. 



2006). However, in other cases the association between the CO and 1.1 mm data is not clear, 
either because of confusion or because the BGPS sources are too compact to identify in ^^CO 
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morphology. Heterodyne follow-up observations using dense gas tracers (NH3, N2H+, HCO^, 
and CS) are being conducted to provide radial velocity, density, chemistry, and temperature 



measurements (Schlingman et al. 2010 Dunham et al. 2010). 



We have detected 8454 clumps in the surveyed fields ( Rosolowsky et al.|2010 ), providing 
a sample of clumps suitable for multi-wavelength and high-resolution studies with existing 
telescopes and future facilities. This sample is free of the usual biases associated with 
observing only regions with signposts of star formation, such as Hll regions or masers. The 
survey depth can be converted into an estimate of limiting mass sensitivity via standard 
estimates, i.e., under the assumption of a temperature, opacity, and distance. Using the 
dust opacity from Ossenkopf & Henning (1994) [k{1.1 mm) = 0.0114 cm^ g~^), a gas-to- 



dust mass ratio of 100, the Bolocam bandcenter from Table |4| and the beam area as given 



in Section 5.1 the mass sensitivity can be written as 



M ^ 



14.3 (1 



,13.0/Ta 



1) 



1 Jy / \ 1 kpc 



D 



Mr, 







(19) 



At the coldest temperature expected (10 K), the mass increases by a factor of 2.9 from the 
mass assuming 20 K. While sensitivity limits will bias the survey against distant, low-mass 
objects, initial results suggest that the Bolocam sources lie at a range of distances, and thus 
range from cores to clumps ( Rosolowsky et ar]|2010 ). We have begun to explore the variation 
of properties with Galactocentric radius with detailed analyses of the BGPS Galactic center 



(Bally et al. 2010) and anti-center (Dunham et al. 2010) data. 



Millimeter-wavelength thermal dust emission reveals the repositories of the densest 
molecular gas, ranging in scale from cores to whole clouds. By pinpointing these regions, 
the BGPS allows the connection of this gas to nascent and ongoing star formation to be 
explored. 



Facilities: CSO (Bolocam) 
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0607793 to the University of Texas at Austin. 
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A. Calculation of Color Corrections 



If an experiment has finite bandwidth (t(i^) 7^ — t'c)), to report a source surface 
brightness at a single frequency, one must assume a source spectrum. The power detected 
from that source is assumed to be 



det 



7] An Io{u)t{u)du 



(Al) 



(Here r] and AQ are the optical efficiency and throughput of the instrument, /o(i^) is the 
nominal (assumed) surface brightness of the source, and t(z/) is the bandpass transmission 
normalized to 1.0 at its peak.) Bolocam is single-moded, meaning that AQ = X^. The 
effective band center z/^ is defined implicitly by the equation 

/ Io{u)t{u)du 



(A2) 



ft{u)du 

In calculating the color corrections and effective band centers, we consider two types of 
source spectra. The most physically motivated is a "greybody" spectrum parametrized as 

iGBiiy) = e{u)B,{T) = (1 - exp (-r,))5.(T) = (1 - exp [-{u / uof])B,{T) (A3) 

where By{T) is the Planck function, e(i/) is the frequency- dependent emissivity, and z/q is 
the frequency at which the optical depth Ty reaches unity. For a frequency much lower than 
the blackbody peak and in the optically thin limit, this reduces to a simple power law 



[V) 



oc 



(A4) 



The color correction K is defined by the relative change in the intensity which would 
be measured at the same effective band center if the source had a different spectrum 
than that assumed for calculating the band center; thus 



J Ii{u)t{u)du 



J lQ{u)t{u)du 



(A5) 



We calculate the effective band centers and color corrections for Bolocam in Table 
4| The fiducial spectrum for quoting the bandcenter (chosen because of its closeness to 



mean Galactic properties; e.g. Reach et al. (1995)) is taken to be that of Equation A3 with 
parameter values T = 20 K, /3 = 1.8, and = 3000 GHz (100 fim). The bandcenters and 
color corrections relative to this fiducial spectrum are also computed for various values of 



a in Equation A4 Over the frequency range of the Bolocam passband, a power law with 
a = 3.5 is a very good approximation to the assumed greybody. 
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B. FITS Header Information 

Below is a sample FITS header from the version 1.0.2 release of the BGPS images 
described in this paper. 



SIMPLE = 
BITPIX = 
NAXIS 
NAXISl = 
NAXIS2 = 
EXTEND = 
ORIGIN = 
IRAF-TLM= 
DATE 

L0NP0LE2= 
LATP0LE2= 
DECFWHM = 
MEANDC = 
STDDC 
UNITS 
PPBEAM = 
CALIB.O = 
CALIB.l = 
CALIB_2 = 
CTYPEl = 
CTYPE2 = 
EQUINOX = 
CDl.l 
CD2_2 
CRPIXl = 
CRPIX2 = 
CRVALl = 
CRVAL2 = 
PV2_1 
PROTITLE= 
CONTACT = 
MAPTYPE = 



T / Written by IDL: Sun Apr 19 17:59:29 2009 
-32 / Bits per pixel 
2 / Number of axes 
1519 / Axis length 
651 / Axis length 

F / File may contain extensions 
NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator 
13:33:03 (13/04/2009)' / Time of last modification 

2009-04-13T19:33:03' / Creation UTC (CCCC-MM-DD) date of FITS header 
180.000000000 /lonpole 
0.00000000000 /latpole 

14.4000 /Deconvolution Kernel FWHM 
2.06882 /Mean DC level 
0.499487 /Std. dev. DC level 
Jy/Beam ' / 

27.7012530153 /pixels per beam 
-3.26472E-15 / 0th coefficient for flux cal (see methods paper 
0.398740 / 1st coefficient for flux cal 
3.32002 / 2nd coefficient for flux cal 
GLON-CAR' / Coordinate Type 

GLAT-CAR' / Coordinate Type 

2000.00 / Equinox of Ref. Coord. 
-0.00199999986216 / Degrees / Pixel 
0.00199999986216 / Degrees / Pixel 

766.378744566 / Reference Pixel in X 
345.530969032 / Reference Pixel in Y 
6.00002238819 / Galactic longitude of reference pixel 
-7.44007626964E-05 / Galactic latitude of reference pixel 
0.00000000000 /Projection parameter 1 
Bolocam Galactic Plane Survey' / 
John Bally' /john.ballyOcolorado . edu 

map ' / 
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BUNIT = 'Jy/Beam 
BMAJ 



0.00988889 / 
0.00988889 / 

/ 



/Units in map 



BMIN 



BPA 



WAVELENG= 



BGPSITER= 



BGPSNPCA= 



BGPSVERS= '1.0 



50 /Iteration number 

13 /number of PGA components subtracted 
/BGPS Processing Version Number 
1.12000 /mm (avoids CO 2-1) 



COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy 
COMMENT and Astrophysics', volume 376, page 359; bibcode 2001A&A. . .376. .359H 

COMMENT Made by the Bolocam Galactic Plane Survey (BGPS) pipeline 

COMMENT described in Aguirre et al 2009 (not yet published) 

COMMENT BGPS data was taken at the Caltech Submillimeter Observatory 

COMMENT Pixel coverage is in the nhitsmap file (each hit represents .Is dwell ti 

COMMENT Pixel weighting is in the weightmap file 

COMMENT Flag counts are in the flagmap file 

COMMENT Deconvolved model is in the model file 

HISTORY PUTAST: Apr 7 09:22:26 2009 World Coordinate System parameters written 

HISTORY Dates and observation numbers included: 

HISTORY 060616_ob3 

HISTORY 060621_o20 

HISTORY 060621_o23 

HISTORY 060622_ol5 

HISTORY 070701_o27 

HISTORY 070702_ol4 

HISTORY 070704_ol6 

HISTORY 070715_ob5 

HISTORY 070706_o22 

HISTORY 070707_o31 

HISTORY 070708_ol5 

HISTORY 070708_ol6 

HISTORY 070709_ol5 

HISTORY 070709_ol6 

WCSDIM = 2 

CDELTl = -0.00199999986216 

CDELT2 = 0.00199999986216 

LTVl = -64. 

LTV2 = -20 . 



LTMl.l = 1 . 

LTM2_2 = 1 . 

WAT0_001= 'system=image' 

WAT1_001= 'wtype=car axtype=glon' 

WAT2_001= 'wtype=car axtype=glat' 



END 
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Fig. 1. — The coverage of the BGPS, showing the continuous coverage in the first quadrant, 
and the regions targeted in the outer Galaxy. The background greyscale is IRAS 100 /xm. 



This preprint was prepared with the AAS I^TjrpC macros v5.2. 
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Fig. 2. — Top: the position and ellipticity of the detectors in the Bolocam focal plane as 
mapped to the sky. The size of the pixel corresponds to its FWHM; it is clear that the 
beams do not completely overlap on the sky. The relative response to source is shown by the 
greyscale, with the darkest being the largest response. The path of the center of each pixel 
as the array is scanned from left to right is shown for a particular orientation of the array 
relative to the scan direction. In this case, the samphng of the sky orthogonal to the scan 
direction is poor. Bottom: the same as the top, but with the array rotated to its optimal 
angle relative to the scan direction. It is clear that the sampling orthogonal to the scan 
direction is now much more uniform. 
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Fig. 3. — The distribution of the pointing cahbration sources across the local sky in Hawaii 
during July 2007 (Epoch V) when the master pointing model was constructed. Note that 
the entire local sky is sampled, and the residual pointing RMS is valid for any point on the 
sky. 
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(a) The pointing model for Epoch V, from which the master reference 
images were derived for subsequent alignment. The left column is 
the pointing model correction and the right the residuals (both in 
arcseconds). The red line indicates the fitted model. 510 pointing 
sources were included, and the final RMS was 5.77" in altitude and 
3.21" in azimuth, or a total RMS offset of 6.6". No systematic offset 
with azimuth was found. 




-10 10 
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(b) The residuals of the pointing model. The Gaussian fits have a = 
5.94, 2.92 in altitude offset and azimuth offset respectively. 



Fig. 4. — The Bolocam pointing model for Epoch V. 
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(a) Left: the Bolocam map in a portion of the I = 30 field. Right: the SLC maps in the same region. 
The green contours on the Bolocam map indicate the positions of the SLC sources. Note that the 
SLC maps only cover a small portion the entire region shown, corresponding to the isolated sources 
indicated by the contours. Negative bowls in both images are forced to zero before cross-correlating. 
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(b) The cross-correlation map of the BGPS and SLC (c) The measured offset between Bolocam and 



maps in Figure 5(a) 



SCUBA Legacy sources (Di Francesco et al. 2008) 



using the method described in Setion 3.4 The circle 
indicates the la region, which is also consistent with 
the error derived from the internal Bolocam pointing 
model. 



Fig. 5. — Comparing the Bolocam positions to those of the SCUBA Legacy Catalog. 
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Fig. 6. — The behavior of the time series as iterative mapping proceeds. 




< 



50 100 150 200 250 300 



Fig. 6a. — Removal of the e and p models. From top top bottom: 1. Raw time series 2. 
After removal of residual 60Hz signal 3. Exponential decay function at scan turnarounds 
subtracted 4. A polynomial (with astrophysical source rejection) is fit to remove large-scale 
variations across the scan. 
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Fig. 6b. — Panel 1 is the same as panel 4 in the previous figure. Before the first cleaning, 
the raw data (black) and remainder data (green) are equal because there is no model. Panel 
2 shows the median-atmosphere subtraction, which is the first- order correction. Panel 3 
shows the astrophysical signal left over from the PCA selection. Panel 4 shows the total 
atmosphere model from PCA and median estimation 
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Fig. 6c. — Same as previous figure, except after 20 iterations. Panel 1 - Tlie deconvolved as- 
trophysical map is returned to a timestream (red) and subtracted from the 'raw' data (black). 
Panel 2 - The median of the remainder (green) is subtracted as the first atmosphere esti- 
mation Panel 3 - The cumulative astrophysical model (red) and the additional astrophysical 
signal from iteration 20 (black) Panel 4 - The atmosphere signal after 20 iterations 
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Fig. 7. — Each PCA component in time can be turned into a map. Above are selected 
images (in inverted greyscale) to illustrate the PCA components and the effect of iteration 
on the process of creating these components. The displays are in Galactic coordinates with 
a color range from -.1 to 1.0 Jy beam~^. The bar in the bottom right indicates a 10' 
length. Top 4 rows: a grid of maps of each of the first 20 PCA components displayed at the 
same scale. It is clear that NCG 7538, the object imaged, has emission in each component. 
There arc varying levels of noise in each component, with the first and second being the 
most obvious atmospheric components. Most of the other noise components are probably 
detector noise correlated among a subset of the bolometers. The streaks arc residuals of the 
scan-turnaround noise that was not removed by the exponential model fit. Bottom 4 rows: 
The same figure, but after 20 iterations. The figure is the breakdown of the atmosphere 
remainder (i.e. the raw data minus the astrophysical model) into eigenf unctions. Very little 
astrophysical flux remains at any level of correlation, though there is some at large (few 
arcminute) scales. 



- 45 - 




Fig. 8. — Left: An illustration of a "glitch" in a single bolometer timestream (red), likely due 
to a cosmic ray strike. Note the acausal ringing due to the application of the downsampling 
filter. The time series for a physically adjacent bolometer is shown in blue, and a bolometer 
which does not pass over the source in green. The PCA atmosphere estimation has already 
been subtracted from the data. Right: The distribution of glitch amplitudes flagged and 
removed from the data in the / = 111 field with their corresponding rate of occurrence. 
Positive-going glitches are the solid histogram and negative-going the dashed. The negative- 
going glitches are likely the ringing from positive-going glitches. Both distributions are 
well-described by a power law above ~ 3 Jy beam~^; below this detection threshold, glitches 
may contribute to the overall noise. 
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Fig. 9. — An illustration of the flagging process using the waterfall plot. Left: A glitch 
is circled in the waterfall plot (bolometer 98, near sample 125). It is most prominent as 
a single hot pixel. Because of the PCA subtraction, however, the effect of the glitch also 
propagates over the other bolometers, appearing as a horizontal stripe; to emphasize this, two 
horizontal dark lines are placed above and below the stripe. Right: The same data displayed 
after flagging out all data affected. Note that the other features in the plot correspond to 
passages of the bolometers over real emission, including NGC 7538, which has a peak of 
~ 7 Jy beam"^. 



-47- 




Fig. 10. — The effect of deconvolution on the iterative process can be seen in its effect on the 
residuals in these images of Sgr B2, all displayed at the same scale. Prom left to right: map, 
residual, model. Top to bottom kernel size: 14'.'4, 21'.'6, 31'.'2, T'.2. The final version of the 
pipeline uses 14'.'4, which has the result of leaving no flux in the residual at the location of 
Sgr B2, and docs not "dig a hole" in the residual map, as the 7'.2 kernel does. The 7'.2 kernel 
also lacks the noise-rejection features of the larger kernels, as can be seen in the bottom right 
panel. 
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Galactic Longitude 

Fig. 11. — The variation of RMS noise as a function of Galatic longitude for the fields in the 
inner Galaxy. The noise reaches a minimum at / = 31 at 11 mJy beam^^, and a maximum 
at / = 57 at 57 mJy beam^^. The shaded region shows the range within which 68% of the 
RMS measurements in a field lie. 
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Fig. 12. — Left: Average calibration curve [V/Jy] versus the mean detector voltage, a proxy 
for atmospheric loading. Asterisks are observations of Mars and crosses of Uranus. The black 
hne is a 2nd-order fit with 0,0 forced (no response if no measurable potential difference) . 
Right: Scaling of the relative response to the atmosphere for a single detector compared 
to the array median. Outliers are from noisy scans that are strongly downweighted. No 
systematic trend in the relative calibration is apparent over the range of atmosphere loading 
observed; the median value is indicated by the dashed grey line. 
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Fig. 13. — An illustration of the relative sensitivity calibration using the atmosphere as a 
calibrator. Black is the median over all bolometers (the Ist-order atmosphere model), red 
and green are individual bolometers before (top) and after (bottom) applying the relative 
calibration. Note the improved agreement. 
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Fig. 14. — Left: The radial profile (solid) of the Bolocam PSF derived from Neptune and 
Uranus observations, including all bolometers. The main lobe is well approximated by a 
Gaussian (dotted) and the sidelobes are due to the Airy pattern from truncation of the 
beam at the cold Lyot stop (dashed). Error bars indicate the azimuthal standard deviation. 
Right: Enclosed flux in an aperture of a given radius based on the measured PSF. This 
curve is the basis of the point source aperture corrections (Table [2]). Vertical dashed lines 
indicated 20, 30, and 40 arcsecond radii. 
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Fig. 15. — The fractional flux recovered as a function of source size for weU-separated 
Gaussian sources with FWHM as indicated. The thin hues show (from top to bottom) the 
effect of increasing the number of PCA components subtracted (3, 7, 10, 13, 16, 21, 26, 
and 31). The thick hne shows the curve for the 13 PCA component cleaning used in the 
released data; the flux recovery drops to 50% at 3'8 for this case. The 3 and 7 PCA cases 
show bumps because atmospheric noise is still present at large scales. The dashed line shows 
the approximation of a a linear high pass fllter with "brick wall" cutoff below 1/10 inverse 
arcminutes. 
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Fig. 16. — Top: the Bolocam 1.1 mm bandpass (thick hne). Also shown 
is the atmospheric transmission at Mauna Kea for 1 mm of precipitable wa- 
ter vapor (thin hne). The grey curve is the MAMBO-2 bandpass (taken 
from 'http : / / www . mpif r-bonn . mpg . de/ div/bolometer/mambo_parameters_020430 . html ) . 
Note that the Bolocam bandpass is well away from the strong water absorption features at 
183 and 325 GHz. Bottom: The effect of line contamination is illustrated with the integrated 
line intensities from the Nummelin et al. (1998) survey of Sgr B2, who found that 22% of 
the flux density in one pointing was due to line emission. Note that the Bolocam passband 
rejects > 90% of the ^^CO flux. However, as shown from the Nummelin et al. (1998) survey, 
SO2 and CH3OH lines can be strong contributors to line flux in the passband due to their 
broad width. Other lines lying in the Bolocam passband include CS(5 — t- 4) and (6 — > 5) 
(245 and 293 GHz), HCN(3 ^ 2) (265 GHz) and HC0+(3 ^ 2) (267 GHz). 
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Fig. 17. — Left: A plot of the values of the pixels in the M07 and BGPS maps after both 
have been convolved to the same resolution on the same grid. The solid line is the expected 
1:1 correspondence; the other lines show the slopes obtained with various lower bounds as 
given in Table [3j There is clearly a great of scatter about the mean. Right: The same 
analysis, but for the M09 SIMBA data. 



Fig. 18. — Examples of the three image types in the IPAC VI. data release. Images are 
zoomed in to show detail: axes are offsets in arcmimites from the center position. Top left: 
MAP. Middle: NOISEMAP Bottom: NHITSMAP. Note the slight "ghosts" in the residual 
(NOISEMAP) at the positions of bright sources. Note also the "basket-weave" pattern in the 
coverage (NHITSMAP) due to gaps in the detector array combined with the scan strategy. 
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Fig. 19. — Images from the BGPS. 
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Fig. 19a. — / = —10.5 to / = 19.5 In this figure and the following, the brightest sources, e.g. 
Sgr B2, Sgr A, and sources near I — 10 and I — 13, appear to be saturated, but this is only 
a display artifact. The astrophysical sources are always much fainter than the atmosphere 
(which is within the dynamic range of the detectors) and therefore do not saturate. The 
noise is more pronounced from I — —1 to I — —2 because this region was observed less. 
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Fig. 19b.— I = 19.5 to I = 49.5. G34.3+0.15, W 51, W 43, W 49, and M 17 appear to be 
saturated, but this is only a display artifact. The 20 < £ < 40 region through the 4-8 kpc 
molecular ring and approximately the termination of the galactic bar is particularly rich in 
clumps. 
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Fig. 19c. — / = 49.5 to / = 74.5. In comparison to the inner galaxy, the 65 < £ < 75 has a 
very sparse population of faint clumps 
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Fig. 19d. — Top: The Cygnus Arm. Note that coverage in 6 is extended to ±1.5 deg. Bottom: 
A zoom-in of the region near DR21, showing at left the data from the M07 IRAM study of 
this region, and at right the BGPS map, in which the filamentary nature of the emission is 



more apparent. Flux densities from the BGPS and M07 are compared in Section 5.5 
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Fig. 19e. — The IC1396 region. In spite of copious CO emission here, there are only two 
faint sources detected in 9 square degrees. 
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Fig. 19f.— Cloud complexes centered at Z = 111 in the Perseus Arm. The NGC 7538 
complex is in the upper-left. 
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Fig. 19g. — W3. The W3(OH)/W3 Main complex is the bright source on the right side of 
the image. 
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Fig. 19h. — W4/5. W5 has one scan performed in RA/Dec instead of Galactic coordinates 
and so has non-uniform noise properties in each square degree 
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Fig. 19i. — Gem OBI. This region has been thoroughly surveyed in NH3, and is discussed 
in Dunham et ah (2010). 
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Table 1. Observing Epochs for the BGPS 



Number 


Begin (UT) 


End (UT) 


I 


2005 Jul 03 


2005 Jul 09 


II 


2005 Sep 05 


2005 Sep 12 


III 


2006 Jun 02 


2006 Jun 30 


IV 


2006 Sep 03 


2006 Sep 19 


V 


2007 Jul 01 


2007 Jul 25 


VI 


2007 Sep 04 


2007 Sep 09 



Table 2. Point Source Aperture Corrections 

RlOAperture Recommended Correction 

40" 1.46 

80" 1.04 

120" 1.00 

object integrated 1.00 



Table 3. Flux comparison with M07 and M09 



Comparison 


Factor 


Factor 


Factor 


Survey 


> 3 MJy sr~^ 


> 10 MJy sr-^ 


> 20 MJy sr-^ 


M07 


1.32 


1.25 


1.21 


M09 


1.51 


1.44 


1.38 



0,4 



0,2 

4-' 

S 0,0 



-0,2 



^1 s. 




-0.4 



3 



0,4 



0.2 r 



0.0 



29,0 
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Fig. 20. — A view of two square degrees of the Galactic Plane centered at / = 28.45, 6 = 0. 
Top: BGPS, Middle: Spttzer-GLIMPSE 8 fim (1'.'9), Bottom: VGPS 20 cm continuum (40"). 
IRDCs appear in the Spitzer negative image as wispy white regions against the grey emission 
background. While some of the Bolocam sources are apparent as IRDCs or strong infrared 
or radio sources, many lack any apparent association with features at these wavelengths. 



The box indicates the zoom-in region in Figures 21 and 22 
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Fig. 21. — A zoom-in of the region shown in Figure 20, showing the VGPS image in inverted 
greyscle, with contours of BGPS emission in black. Selected sources from the RIO catalog 
are marked with X's and their catalog numbers. Note that while Hll regions appear at or 
near the peaks of millimeter emission for clumps 3913, 3922, 3935, 3955, and 3971, others 
have no associated emission (3897, 3921, 3923). 
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Fig. 22. — A zoom-in of the region shown in Figure 20 , showing the Spitzer-GLlMPSFj image 
in inverted greyscle, with contours of BGPS emission in black. IRDCs appear as white in 
this image. Selected sources from the RIO catalog are marked with X's and their catalog 
numbers. Note the high degree of correspondence between the IRDC morphology and the 
millimeter emission for BGPS sources 3923 and 3955. 
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Fig. 23. — The filling factor of millimeter continuum emission from the BGPS compared to 
C0(1 — i- 0) emission from Jackson et al. (2006). The axis on the left side gives the filling 
fraction above 3a of the local RMS from the BGPS (solid line). The right axis gives the CO 
filling fraction above 50 (dotted) and 100 (dashed) K km s~^. 
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Table 4. Comparison of band centers, color corrections, and flux ratios for Bolocam and 

MAMBO. 





Bolocam 


K 




MAMBO 


'S'bolocam/ 




Vc [GHz] 


Eq. 


A5 


Vc [GHz] 


5'mambo^ 


c 


271.1 


1.000 


262.0 


1.124 


1.0 


267.8 


0.98 


8 


254.8 


1.051 


1.5 


268.6 


0.98 


6 


256.4 


1.072 


2.0 


269.4 


0.98 


7 


258.0 


1.090 


2.5 


270.0 


0.990 


259.5 


1.104 


3.0 


270.7 


0.996 


261.1 


1.115 


3.5 


271.4 


1.004 


262.5 


1.124 


4.0 


272.2 


1.015 


264.0 


1.130 


4.5 


272.9 


1.030 


265.4 


1.134 


5.0 


273.8 


1.049 


266.8 


1.137 



^The spectral index as deflned by Equation A4 



b 



The expected ratio of flux densities, assuming 



identical apertures, the passbands of Figure [I6| and 
the SED indicated. 

^The flducial SED is taken to be that of Equation 



A3 with T = 20 K, /3 = 1.8, and unity optical depth 
at 100 /im. 



